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’9  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

A  tunnel  diversion  project  has  been  authorized  that  will  divert  the  full  flow  of  the 
Clover  Fork  near  Harlan,  KY,  to  an  outlet  downstream  of  the  city.  The  project  has  a 
diversion  structure  at  the  upstream  end  of  the  project  that  will  direct  the  flow  of  the 
Clover  Fork  into  four  inverted  U-shaped  tunnels  that  average  1,936  ft  in  length.  The  narrow 
floodplain  upstream  of  the  project  site  is  filled  with  debris  that  could  potentially  reach 
the  entrance  to  the  tunnels  during  the  Standard  Project  Flood  (SPF) . 

A  1: 30-scale  physical  model  was  constructed  to  assess  the  tunnel  entrance  for  the 
improvement  of  hydrodynamic  conditions  and  the  mitigation  of  potential  floating  debris 
problems.  A  two-dimensional  depth-averaged  numerical  hydrodynamic  model  was  used  in 
conjunction  with  the  physical  model  to  screen  alternatives  that  would  streamline  the 
approach.  The  physical  model  investigated  these  alternatives  for  their  capability  to  pass 
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debris  by  reproducing  the  SPF  hydrograph  while  simultaneously  introducing  predicted  volume 
of  scaled  debris . 

Features  reproduced  in  the  physical  model  included  about  2,000  ft  of  the  natural 
channel  upstream  of  the  tunnel  entrance,  the  diversion  structure,  the  proposed  tunnel 
entrance,  and  450  ft  of  the  tunnels.  The  testing  resulted  in  two  entrance  designs  that  not 
only  passed  the  SPF  flow  with  debris  but  also  provided  sufficient  freeboard  along  the 
diversion  structure. 


PREFACE 


The  model  investigation  and  analysis  described  herein  were  authorized  by 
Headquarters,  US  Army  Corps  of  Engineers  (USACE) ,  on  24  September  1987  at  the 
request  of  the  US  Army  Engineer  District,  Nashville  (ORN) .  The  model  tests 
were  conducted  on  an  accelerated  schedule  at  the  request  of  ORN  during  the 
period  May  1988  to  December  1988  in  the  Hydraulics  Laboratory  of  the  US  Army 
Engineer  Waterways  Experiment  Station  (WES) . 

Testing  was  conducted  under  the  direction  of  Messrs.  F.  A.  Herrmann,  Jr., 
Chief  of  the  Hydraulics  Laboratory;  R.  A.  Sager,  Assistant  Chief  of  the 
Hydraulics  Laboratory;  G.  A.  Pickering,  Chief  of  Hydraulics  Structures 
Division  (HSD) ;  and  J.  F.  George,  Chief  of  Locks  and  Conduits  Branch  (LCB) , 

HSD.  Physical  model  tests  were  conducted  by  Ms.  Sandra  K.  Martin  and  Mr.  Van 
Stewart,  LCB.  The  computer  code,  STREMR,  used  in  this  study  was  written  by 
Dr.  Robert  S.  Bernard,  Reservoir  Water  Quality  Branch,  HSD.  Dr.  Bernard's 
expertise  and  timely  completion  of  the  model  tests  were  an  integral  part  of 
the  successful  completion  of  this  study.  This  report  was  prepared  by 
Ms.  Martin  and  edited  by  Ms.  Marsha  C.  Gay,  Information  Technology  Laboratory, 
WES. 

Prior  to  model  testing,  Ms.  Martin  visited  the  proposed  project  site 
with  ORN  personnel  in  Harlan,  KY,  and  vicinity  to  inventory  debris  in  the 
basin.  Visitors  to  WES  during  the  model  tests  included  Messrs.  Sam  Powell  and 
Tom  Munsey,  USACE;  Mr.  Lynn  Richardson,  US  Army  Engineer  Division,  Ohio  River; 
LTC  William  Allen,  Deputy  District  Engineer,  ORN;  Messrs.  Warren  Bennett, 

Euclid  Moore,  Hank  Phillips,  Tom  Allen,  Jody  Stanton,  Ben  Couch,  and  Marvin 
Simmons,  ORN,  and  several  representatives  from  private  industry. 

Acting  Commander  and  Director  of  WES  during  preparation  of  this  report 
was  LTC  Jack  R.  Stephens,  EN.  Technical  Director  was  Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC)  UNITS 
OF  MEASUREMENT 


Non-SI  units  of  measurement 
(metric)  units  as  follows: 

_ Multiply _ 

cubic  feet 
degrees  (angle) 
feet 

miles  (US  statute) 
square  feet 
square  miles 


used  in  this  report  can  be 


_ Bx _ 

0.4535924 

0.0 1745329 

0.3048 

1.609344 

0.09290304 

1.609344 


converted  to  SI 

To  Obtain _ 

kilograms 

radians 

metres 

kilometres 

square  metres 

square  kilometres 
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Figure  i.  Location  and  vicinity  map 
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CLOVER  FORK  TUNNEL  DIVERSION  PROJECT 
HARLAN.  KENTUCKY 
Hydraulic  Model  Investigation 

PART  I:  INTRODUCTION 

The  Prototype 

1.  Martins  Fork  joins  Clover  Fork  in  Harlan,  KY,  approximately  1  mile* 
south  of  Clover  Fork's  confluence  with  Poor  Fork  and  the  Cumberland  River 
(Figure  1) .  Historical  flooding  of  the  three  forks  can  be  characterized  as 
fast  rising  (8-10  hours)  and  short  lived  in  duration  (approximately 
24-30  hours).  Both  the  Martins  Fork  and  Clover  Fork  channels  are  approxi¬ 
mately  100  ft  wide,  have  stream  slopes  of  about  11  ft  per  mile,  and  have  banks 
5  to  15  ft  high.  Clover  Fork  and  Martins  Fork  drain  approximately  104  and  117 
square  miles,  respectively,  of  primarily  forested  land.  Clover  Fork's  narrow 
floodplain,  averaging  approximately  1,000  ft  in  width,  is  cluttered  with 
trash,  dead  brush,  potentially  unstable  housing  such  as  trailers,  and 
stockpiles  of  logs  from  clearing  operations  in  the  area.  The  shallow  channel 
banks  show  signs  of  instability  and  are  lined  with  large  trees  that  lean 
precariously  toward  the  channel  bed. 

2.  A  tunnel  diversion  project  is  authorized  that  will  provide  Standard 
Project  Flood  (SPF)  protection  by  diverting  the  full  flow  of  Clover  Fork 
around  the  Harlan  central  business  district.  A  diversion  structure  will  be 
constructed  at  the  upstream  end  of  the  project  area  to  direct  the  flow  from 
Clover  Fork  into  four  tunnels  that  discharge  downstream  of  Harlan  (Plate  1) . 
The  size  of  the  Harlan  tunnels  and  the  diversion  levee  height  were  designed  to 
ensure  SPF  protection  with  an  appropriate  level  of  freeboard  on  the  diversion 
structure.  Each  tunnel  has  an  inverted  U-shaped  cross  section,  approximately 
33  ft  wide  by  32  ft  high.  Separation  between  tunnels  is  approximately  35  ft. 
The  tunnels,  averaging  1,936  ft  in  length,  are  designed  to  pass  a  maximum  SPF 
discharge  of  52,800  cfs  with  an  assumed  percentage  of  debris  blockage  at  the 
entrance.  The  proposed  headwater  elevation  that  determined  the  levee  height 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  3. 
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was  ilculated  assuming  an  entrance  loss  coefficient  due  to  blockage  by 
debris . 


Purpose  and  Scope  of  Model  Investigation 

3.  One  of  the  major  project  concerns  was  the  potential  retention  of 
floating  debris  at  the  upstream  face  of  the  tunnels.  Blockage  at  the  portal 
face  resulting  from  the  accumulation  of  debris  would  raise  the  headwater  in 
the  approach  area  and  thereby  decrease  the  available  freeboard  on  the  diver¬ 
sion  structure.  The  study  was  initiated  to  evaluate  the  adequacy  of  the 
proposed  tunnel  entrance  to  pass  debris  and  to  maintain  the  design  freeboard, 
and/or  to  develop  modifications  that  would  not  only  improve  the  hydrodynamic 
conditions,  but  also  mitigate  potential  debris  problems  at  the  upstream 
portal . 

4.  Both  a  physical  and  a  mathematical  model  were  used  to  evaluate  the 
proposed  entrance  approach.  The  mathematical  model  was  used  in  conjunction 
with  the  physical  model  to  screen  alternatives,  while  the  physical  model 
investigated  these  alternatives  in  detail  for  their  ability  to  pass  debris  by 
simultaneously  reproducing  the  SPF  hydrograph  while  an  expected  volume  of 
scaled  debris  was  introduced.  The  main  objectives  of  the  study  were  to 
evaluate  and  improve  the  hydrodynamic  design  of  the  entrance  approach,  and  to 
select  a  design  configuration,  especially  regarding  the  tunnel  piers  and 
portal  face,  that  minimizes  head  losses  associated  with  debris. 
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PART  II:  THE  MODELS 


Physical  Model 


5.  The  1: 30-scale  physical  model  reproduced  approximately  2,000  ft  of 
the  natural  channel  upstream  of  the  entrance,  the  diversion  structure,  the 
proposed  channel  approach,  and  450  ft  of  the  tunnels  (Figure  2).  All  elements 
of  the  model  except  the  tunnels  were  constructed  by  molding  sand  and  cement 
mortar  to  sheet  metal  templates.  The  tunnels  were  constructed  of  transparent 
plastic  for  observation  of  water  and  debris  flow  conditions. 


Mode  1  nppur t enances 


6.  Water  used  in  the  operation  of  the  model  was  supplied  by  a  circulat¬ 
ing  .system.  Discharges  were  measured  with  a  venturi  meter  installed  in  the 
flow  line  and  were  baffled  before  entering  the  model.  Steel  rails  graded  to 
specific  elevations  were  placed  alorg  both  sides  of  the  model  to  serve  as 
supports  for  measuring  devices.  Velocities  were  measured  with  an  electronic 
velocity  meter  mounted  to  permit  mea  urement  of  flow  from  any  direction  and  at 
any  depth.  Water- surf ace  elevations  were  measured  with  point  gages,  and 
pressures  were  measured  with  piezometers.  Tailwater  conditions  were  regulated 
by  adjusting  gates  located  at  the  downstream  end  of  each  tunnel  and  monitoring 
piezometers  located  on  the  floor  of  the  tunnels.  Different  designs  for  SPF 
peak  flow  conditions  and  SPF  hydrographs  with  debris  were  recorded 
photographically . 

Scale  relations 

7.  The  equations  of  hydraulic  similitude,  based  on  the  Froude  criteria, 
were  used  to  express  mathematical  relations  between  the  dimensions  and 
hydraulic  quantities  of  the  model  and  prototype.  General  relations  for 
transferring  model  data  to  prototype  equivalents  are  as  follows: 


wharac  t  e  r  l  s r.  ic 


Dimension* 


Scale  Relations 
Mode  1  ,  Procot  vi>e 
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Looking  upstream  b.  Looking  downstrea: 


Sl'REMR,  the  numerical  model  selected  tor  this  studv.  is  a  finite- 
diltii'ince  computer  code  that  uses  the  governing,  equations  tor  depth-averaged 
i  nci'inpi't  ss  i  hi  t  1 1  ow  in  two  dimensions  (d-D).  The  model,  was  developed  bv 
L'r.  Robert  S.  Bernard  tor  the  Repair,  Evaluation.  Maintenance  and 
Rehabi 1 i t  at  ion  < REMR  )  Research  Program  sponsored  bv  Headquarters.  US  Armv 
I'orps  ot  Engineers  (USAGE)  (Bernard  148/ ).  Lateral  boundaries  and  obstacles 
c.iti  have  anv  shape,  and  bathvmetrv  can  varv  with  position.  The  code  employs  a 
predicto^ -corrector  scheme  to  solve  the  discretized  governing  equations; 
resistance  on  the  bottom  is  computed  from  Manning's  equation.  The  tlow  mav  be 
either  steadv  or  t ime -dependent . 

u .  The  fundamental  assumption  in  the  STREMR  code  is  that  the  Mow 
upstream  ot  a  hvdraulic  structure  can  be  described  bv  the  depth-averaged 
equations  ot  motion  tor  an  incompressible  tluid.  Additionally,  it  is  assumed 
that  tlte  displacement  ot  the  free  surface  from  its  initial  elevation  is  small: 
therefore-,  the  computational  free  surface  becomes  a  rigid  lid.  This  assump¬ 
tion  is  exact  at  a  Froude  number  equal  to  U.  but  gradually  breaks  down  as  the 
1'roude  number  approaches  L. 


PART  III: 


TESTS  AND  RESULTS 


10.  Seven  different  approach  designs  were  constructed  in  the  physical 
model.  To  streamline  the  approach  hydrodynamically ,  evaluate  passage  of 
debris,  and  assess  entrance  losses  at  the  portal  face,  three  kinds  of  tests 
were  conducted.  In  the  first  kind  of  test,  velocities  and  water  -  surface 
profiles  were  used  to  compare  types  1-7  design  approaches  during  the  SPF  peak 
flow.  ^econdl",  the  response  of  types  4-7  approaches  to  floating  debris  was 
tested  by  reproducing  the  SPF  hydrograph  while  introducing  debris.  Finally, 
the  coefficient  for  entrance  loss  was  determined  for  variable  blockage 
conditions  on  the  type  6  design  approach. 

Approach  Channel  Designs 


11.  The  seven  different  approach  designs  are  described  as  follows: 

a.  Type  1  design  approach.  Original  design  approach  channel  with 
the  portal  face  having  a  vertical  slope  of  4V : 1H  without  tunnel 
piers  (Photo  1,  Plate  2). 

b.  Type  2  design  approach.  The  same  approach  layout  as  the  type  1 
design  approach  except  that  the  portal  face  was  vertical,  a 
5-ft-radius  bellmouth  was  installed  on  the  tunnel  openings,  and 
three  tunnel  piers  and  two  bridge  piers  were  installed  accord¬ 
ing  to  drawings  provided  by  the  US  Army  Engineer  District 
(USAED) ,  Nashville  (Photo  2,  Plate  3). 

c.  Type  3  design  approach.  Same  as  the  type  2  design  approach 
except  that  the  existing  right  bank  approach  wall  was  replaced 
with  a  curvilinear  transition  section  (Photo  3,  Plate  4). 

d.  Type  4  design  approach.  The  type  3  design  approach  was  modi¬ 
fied  by  lowering  the  noses  on  tunnel  piers  1  and  3,  remov¬ 
ing  the  bridge  piers  on  top  of  tunnel  piers  1  and  3,  and 
replacing  them  with  one  bridge  pier  located  on  the  center  of 
tunnel  pier  2  (Photo  4,  Plate  5). 

e.  Type  5  design  approach.  The  tops  of  the  tunnel  piers  in  the 
type  4  design  approach  at  the  portal  face  were  raised  to  ele¬ 
vation  1210.0,*  and  the  sloping  pier  noses  were  extended  up 
to  this  elevation  (Photo  3,  Plate  6). 

f.  Type  6  desip.n  approach.  The  type  5  tunnel  piers  were  com¬ 
pletely  removed  and  replaced  with  piers  designed  similar 


*  All  elevations  (el)  cited  herein  are  in  feet  referred  to  the  National 
Geodetic  Vertical  Datum  (NGVD) . 
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to  those  in  Engineer  Manual  (EM)  1110-2-1601  (Headquarters, 
USACE,  1970)  (Photo  6,  Plate  7).  The  5-ft-radius  bellmouth  on 
each  tunnel  was  replaced  with  a  10. 5- ft  bell. 

g.  Type  7  design  approach.  The  final  pier  design,  furnished  by 
Nashville  District,  consisted  of  120-ft-long  tunnel  piers  with 
a  sloping  nose  that  had  an  8-ft  radius  (Photo  7,  Plate  8). 

Since  the  top  of  the  piers  extended  above  the  headwater  ele¬ 
vation  along  the  roadway  alignment,  the  bridge  pier  was  not 
necessary . 

In  types  2-7  designs,  a  warped  transition  section  was  needed  to  tie  the  tunnel 
piers  and  the  bellmoulli  smoothly  into  the  tunnel  walls. 

12.  The  bridge  and  tunnel  pier  designs  in  types  2,  3,  4,  and  5  were 
provided  by  the  Nashville  District.  The  bellmouth  entrance  to  the  tunnels, 
the  curvilinear  right  bank  wall,  and  the  tunnel  pier  design  in  type  6  were 
recommended  for  testing  as  modifications  for  the  improvement  of  hydraulic 
conditions.  The  bridge  piers  were  needed  to  support  a  roadway  proposed  to 
cross  the  tunnel  piers  at  an  alignment  specified  in  the  Harlan  Diversion 
Project  Feature  Design  Memorandum  (USAED,  Nashville,  1988).  Type  7  tunnel 
piers  were  selected  for  testing  from  three  options  designed  by  the  staff  at 
Nashville  District  that  potentially  satisfied  both  hydraulic  and  geotechnical 
attributes . 


Model  Data  and  Observations 


13.  Initially,  there  was  some  concern  that  the  abrupt  transition  from 
the  natural  channel  into  the  tunnels  would  skew  the  flow  distribution  down 
each  tunnel,  thereby  decreasing  the  total  capacity  of  the  tunnels.  The  tun¬ 
nels  were  sized  assuming  an  even  distribution  of  flow  down  each  tunnel  and  a 
level  headwater  ^ross  the  portal  face.  Surface  currents,  as  exhibited  by 
confetti,  displayed  a  tendency  to  pull  to  the  left  bank  wall.  Therefore, 
velocities  were  taken  to  ensure  that  placement  and  shape  of  piers  did  not 
drastically  alter  the  distribution  of  the  flow  down  each  tunnel.  Water- 
surface  elevations  at  the  portal  face  were  taken  to  evaluate  the  hydraulic 
efficiency  of  the  various  entrance  conditions.  Results  of  these  tests  are 
summarized  in  Table  1. 

14.  Velocities  and  water-surface  profiles  were  recorded  for  types  1-7 
design  approaches  during  steady-state  flow  conditions.  The  peak  inflow  of 
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52,800  cfs  was  reproduced  with  a  corresponding  peaK  tailwater  elevation  of 
1195.6  at  the  downstream  tunnel  outlet.  The  tailwater  was  reproduced  by 
controlling  gates  at  the  outlet  to  each  tunnel  in  the  model  until  piezometers 
located  420  ft  (prototype)  downstream  from  the  tunnel  entrance  produced  a 
point  on  the  hydraulic  grade  line  equivalent  to  el  1198.7  at  that  station. 

The  velocities  were  taken  at  five  elevations  for  stations  2.5  ft,  45  ft,  and 
90  ft  upstream  of  the  portal  face  and  along  the  center  line  of  the  tunnels. 

At  all  other  stations,  velocities  were  taken  at  near-bottom,  near- top,  and 
middepth  elevations  laterally  across  the  section.  Photos  8  and  9  reflect 


current  patterns  for  type  1,  the  original  design,  and  type  7,  the  selected 
design,  respectively.  These  photographs  also  indicate  the  close  correlation 
between  STREMR  results  discussed  in  the  following  section  and  physical  model 
observations.  Velocities  measured  with  the  type  1  design  approach  and  type  7 
design  approach  are  shown  in  Plates  9  and  10,  respectively. 

Numerical  Model  Tests 

15.  The  STREMR  2-D  model  was  used  simultaneously  with  the  physical 
model  to  assess  the  entrance  conditions  to  the  proposed  tunnel  site.  The 
model  grid  was  originally  developed  for  a  far-field  approach  to  assess  general 
flow  patterns  in  the  natural  channel  for  approximately  3,000  ft  upstream  of 
the  tunnel  entrance.  Figure  3  shows  the  general  layout  grid,  streamlines,  and 
velocity  vectors. 

16.  The  grid  was  then  refined  for  an  intermediate  range  covering 
approximately  1,000  ft  upstream  of  the  portal  face.  Model  tests  performed 
with  this  grid  verified  the  location  of  local  eddies  in  the  physical  model, 
identified  flow  separation  points,  and  approximated  the  direction  and  relative 
magnitude  of  the  velocity  vectors  observed  in  the  physical  model.  As  an 
example  of  the  correlation  between  models,  Figure  4  displays  velocity  vectors 
measured  in  the  physical  model  and  those  calculated  from  an  intermediate- range 
grid  of  the  type  2  design  approach.  The  type  3  design  approach,  in  which  the 
right  bank  approach  wall  was  modified,  was  recommended  as  a  result  of  these 
tests  (Plates  11-13). 

17.  Finally,  two  grids  were  developed  for  close -approach  calculations 
with  the  type  4  design  approach  piers  in  place:  one  with  the  original  right 
bank  approach  and  one  with  the  type  3  curvilinear  right  bank  approach.  As  can 
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be  seen  bv  the  streamlines  and  velocity  vectors  shown  in  Plates  14-16, 
smoothing  the  right  bank  wall  removed  the  separation  point  and  eddy  observed 
bv  the  original  right  bank  wall,  In  the  physical  model,  the  flow  separation 
was  eliminated  as  was  the  eddy;  however,  some  local  turbulence  occurred  at  the 
bank-water  interface.  While  the  correlation  between  the  models  was  good, 
three-dimensional  effects  cannot  be  reproduced  in  a  2-D  model.  The  mathemati¬ 
cal  model  also  supported  conclusions  drawn  from  the  velocity  data  in  the 
physical  model,  that  each  tunnel  carries  approximately  one-quarter  of  the 
total  discharge . 


Debris  Tests 


Vo  1 ume 

18.  Since  the  major  concern  of  the  model  study  was  to  assess  the 
ability  of  the  various  designs  to  pass  debris  through  the  tunnels  and  max¬ 
imize  available  freeboard  on  the  diversion  structure,  a  testing  methodology 
was  developed  in  which  a  specific  volume  of  scaled  debris  was  introduced  into 
the  model  during  the  SPF  hydrograph.  The  volume  of  debris  expected  to  reach 
the  upstream  entrance  to  the  tunnels  during  a  storm  of  this  magnitude  was 
determined  by  Nashville  District.  The  estimated  volume  of  debris  was 
1,923,000  ft3,  or  approximately  72  ft3  in  the  model.  This  estimate  was  based 
on  information  obtained  from  a  field  investigation  in  April  1988,  available 
literature  concerning  inventory  methodologies  (McFadden  and  Stallion  1976), 
and  photographs  taken  during  historical  flood  events  in  the  region. 

Mixture 

19.  The  debris  accessible  to  transport  in  the  Clover  Fork  floodplain  is 
a  heterogeneous  mixture  having  various  shapes,  specific  weights,  and  sizes. 

For  purposes  of  quantifying  the  debris,  it  was  broken  into  several  categories: 
large  trees,  small  trees,  branches  and  twigs,  housing,  and  housing  parts.  A 
worst-case  distribution  by  volume  was  to  consist  of  70  percent  large  objects 
(defined  as  those  with  nominal  dimensions  exceeding  35  ft)  including  intact 
houses  and  trailers.  The  remaining  30  percent  was  to  include  smaller  objects 
such  as  broken  trailers,  twigs,  and  branches.  A  sensitivity  test  using  this 
debris  scenario  on  the  type  4  design  approach  during  the  SPF  hydrograph 
resulted  in  overtopping  of  the  diversion  structure.  The  elevation  of  the 
structure  in  the  model  is  1210.0.  The  final  design  top  elevation  of  the 
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diversion  structure  is  approximately  1212.0.  Skewing  the  distribution  of 
debris  to  the  larger  objects  was  felt  to  be  inappropriate  by  personnel  from 
USAGE,  the  US  Army  Engineer  Division,  Ohio  River,  Nashville  District,  and  the 
US  Army  Engineer  Waterways  Experiment  Station.  A  more  reasonable  distribution 
was  agreed  upon  in  which  80  percent,  of  the  total  volume  consisted  of  smaller 
objects  including  broken  pieces  of  trailers  and  housing  of  variable  buoyancy, 
as  well  as  twigs  and  branches.  The  remaining  20  percent  included  trees  from 
35  to  90  ft  tall,  some  of  which  were  fabricated  to  have  nonbuoyant  root  sys¬ 
tems.  Table  2  lists  the  distribution  of  the  debris. 

The  hvdrograph  with  debris 

20.  The  inflowing  SPF  hydrograph  used  during  all  tests  began  3  hours 
prior  to  the  peak  and  continued  through  the  peak  for  3  hours  on  the  recession 
limb  for  a  total  of  approximately  6  hours,  prototype.  Consequently,  the  be¬ 
ginning  flow  was  47,000  cfs,  rising  to  a  peak  of  52,650  cfs,  and  receding  to 
43,500  cfs.  Simultaneously,  the  tailwater  hydrograph  was  reproduced  by  opera¬ 
tion  of  the  downstream  gates  on  the  tunnels.  Both  the  inflowing  hydrograph 
and  the  tailwater  hydrograph  were  provided  by  Nashville  District  (Figure  5). 
The  debris  was  introduced  evenly  throughout  the  hydrograph  except 


Figure  5.  SPF  discharge  and  stage  hydrographs 
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HYDRAULIC  GRADELINE  ELEVATION 


during  two  instances  when  upstream  "bridge  failures"  were  simulated  by  dumping 
a  specific  volume  (137,500  ft  )  of  debris  into  the  model  at  once.  These 
failures  occurred  approximately  1  hour  prior  to  the  peak  and  during  peak 
discharge.  Ivhile  numerous  scenarios  of  both  debris  inflow  and  the  inflowing 
hydrograph  could  have  been  tested,  the  specific  conditions  selected  were 
assumed  to  be  representative  and  tested  on  types  4-7  design  approaches  to 
evaluate  the  relative  differences  between  designs.  Debris  tests  were  not 
conducted  on  design  types  1-3. 

Test  sequence 

21.  The  SPF  hydrograph,  tailwater  hydrograph,  and  debris  conditions 
described  in  the  preceding  paragraph  were  run  with  types  4,  5,  6,  and  7  design 
approaches.  Design  performance  was  evaluated  by  observation  of  debris  accumu¬ 
lation  at  the  tunnel  entrance  and  the  maximum  water- surface  elevation  in  the 
approach  area.  Time-lapse  photographs  were  taken  at  5-min  intervals  in  the 
model  and  a  video  camera  recorded  the  entire  debris  test.  Photos  10-13  taken 
at  approximately  2  hours  past  the  peak  SPF  discharge  during  each  test  show  the 
relative  difference  in  types  4-7  designs,  respectively.  The  elevation  lines 
painted  on  the  model  along  the  diversion  structure  were  used  to  observe  the 
maximum  headwater  during  the  test.  The  peak  headwater  elevation  and  total 
losses  for  each  test  condition  are  found  in  Table  2.  Further  documentation  of 
the  debris  tests  can  be  found  in  Martin  (1989). 

Hydraulic  Losses 


Theory  and  equations 

22.  A  portion  of  the  testing  program  was  devoted  to  the  evaluation  of 
the  losses  at  the  tunnel  entrance.  Prior  to  the  model  tests,  the  tunnels  were 
sized  assuming  an  entrance  loss  coefficient  Kg  of  1.02.  This  coefficient 
was  based  on  a  30  percent  blockage  and  determined  according  to  the  following 
equation : 


whe  re 


K„  =  0.5 
E 


/  at 

\%0PEN  x  A 


o 

A-j,  =  tunnel  area  “  922.7  ft 
%0PEN  =  70  percent  (30  percent  blockage  assumed) 


(1) 
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The  sum  of  the  Loss  coefficients  K^,  is  applied  to  the  velocity  head  at  the 
exit  to  determine  the  total  head  loss  H-p  from  downstream  to  upstream  for 
tunnels  flowing  full  with  outlet  control  as  follows: 


2 

Hx  -  (Kg  +  S  +  V  x  ^  (2) 

where 

Kp  =  loss  due  to  friction  =  (29.1  n^D/R^7^  where  n  =  Manning's 

coefficient.  0.018;  L  =  average  length  of  tunnels,  2.100  ft;  and 
R  =  hvdraulic  radius.  8.03  ft 

Kq  =  outlet  loss  coefficient  =  1  (assumed) 

v  =  exit  velocity  =  Q/A  =  52,800/(4  x  922.7)  =  14.3  fps 

g  =  gravitational  constant  =  32.16  ft/sec ^ 

and 

Kt  =  K£  +  Kp  +  KQ  (3) 

During  assumed  design  conditions,  when  Kp  equals  1.23  and  Kp  equals  1.02. 
Hj  is  equal  to  10.3  ft.  Adding  10.3  ft  to  the  design  tailwater  elevation  of 
1195.6  gives  a  headwater  elevation  of  1205.9.  While  the  tunnels  were  sized 
based  on  30  percent  blockage,  the  elevation  of  the  diversion  structure  was 
raised  to  accommodate  headwater  conditions  for  a  50  percent  blockage  factor. 

23.  Equation  2  assumes  no  velocity  head  in  the  upstream  approach  when 
in  actuality  a  velocity  head  does  exist.  To  calculate  the  entrance  loss 
coefficient  Kp  for  the  various  headwater  conditions  observed  during  testing, 
it  becomes  necessary  to  incorporate  the  upstream  velocity  head  into  Equation  2 
using  the  energy  equation, 

2  2 

v  v 

1  2 

—  +  HW  =  TW  +  —  +  losses  (4) 

2g  2g 


where 

Vp  =  6.5  fps  (average  velocity  upstream  of  the  tunnels,  taken  from  model 
data ) 

HW  =  headwater  elevation  observed  during  testing 
TW  =  tailwater  =  1195.6 

-  exit  velocity  =  14.3  fps 
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If  v~  2g  +  losses  is  replaced  with  the  right  side  of  Equation  2.  the  follow¬ 
ing  is  obtained: 

2  2 

V1  V2 

+  (K£  +  Kf.  +  Kq)  X  —  (b) 

Defining  as  HW  minus  TV  .  using  the  friction  coefficient  Kp  =  1.18  for 

the  model,  and  substituting  the  values  assumed  for  vp  ■  v2  '  an<^  ' 

Equation  o  can  be  solved  tor  Kp  .  (The  prototype  n  value  of  the  plastic  in 
the  model  is  0.016.  which  is  slightly  lower  than  the  assumed  n  value  of 
shotcrete,  0.018.  resulting  in  a  lower  coefficient  of  friction  for  use  in  the 
model  equat ion .  ) 

KE  -  0.3145  Ht  -  1.973  (6) 


Losses  due  to  debris 

24.  Table  3  presents  the  loss  coefficients  Kp  and  K^,  for  observed 
headwater  conditions  in  the  model  during  debris  testing.  For  comparative 
purposes  the  table  also  contains  the  estimated  (EST)  design  conditions  based 
on  Equations  1  and  2.  For  all  designs,  except  tvpes  4  and  5  with  debris,  the 
actual  conditions  result  in  lower  headwater  elevations  than  those  estimated  by 
the  design  criteria.  Tvpes  6  and  7  resulted  in  the  lowest  headwater  elevation 
for  conditions  with  and  without  debris.  A  chart  for  comparing  head  loss 
measurements  is  provided  in  Figure  6. 

Losses  due  to  blockage 

25.  Another  set  of  tests  was  conducted  to  determine  the  head  loss 
associated  with  various  blockage  conditions.  Five  tests  were  conducted  bv 
blocking  the  upper  portion  of  the  tunnels  in  the  type  6  design  approach  with  a 
solid  nonporous  cover  set  at  elevations  that  blocked,  respectively,  10,  20. 

30,  40,  and  50  percent  of  the  cross-sectional  area  of  each  tunnel.  The 
results  are  as  follows: 


Blockage 
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Headwater 

El 

ht  .  ft 

KE 

kt 

10 

1202.6 

7.0 

0.23 

2.41 

20 
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0.35 

2.53 

30 

1204.5 

8.8 

0.80 

2.98 

40 

1208.9 

13.2 

2.18 

4.36 

50 

1217.0 

21.4 

4.76 

6.94 
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[^2!]  WITH  DEBRIS  [\2S]  NO  debris 
Figure  6.  Comparison  ot  head  loss  conditions 

26.  Additional lv .  six  rests  were  conducted  bv  covering  the  intakes  with 
variable  porositv  stainless  steel  screens.  These  screens  had  wire  diameters  ana 
openings  which  produced,  respec t ivei v .  23.6.  29.4.  34.4,  39.8.  44.6.  and  51.1 
percent  areal  blockages  of  the  tunnel  entrances.  The  results  of  these  tests  are 
as  follows: 


Blockage  Headwater 


Percent 

El 

.  1 1 

ke 

kt 

2  3.6 

1202.8 

/  .2 

0.29 

2.4/ 

29.4 

1202.8 

7.2 

0.29 

2.4/ 

34.4 

1203.8 

8.2 

0.61 

2  .  /  9 

CO 

CD 

oc 

1204.0 

8.4 

0.6/ 

2.85 

44 . 6 

1204.8 

9.2 

0.92 

3  .  10 

51  .  1 

1209.0 

13.4 

2.24 

4.42 

The  plots  in  Figure  /  show  percent  blockage  versus  head  loss  for  porous  and 
nonporous  materials. 
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BLOCKAGE,  PERCENT 


Figure  7.  Comparison  of  head  loss  for  porous  and  nenporous 
blockages  using  the  type  6  design  approach 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


2  7 .  Results  of  the  physical  and  numerical  model  studies  led  to  a  number 
of  conclusions.  Where  possible,  all  transition  sections  from  the  natural 
channel  to  the  entrance  channel  that  are  at  or  below  the  design  headwater 
elevation  should  be  curved  to  prevent  flow  separation  and  local  eddies  from 
developing  (for  example,  the  right  wall  transition  in  the  type  3  design 
approach) . 

18.  The  3-ft-radius  belimouth  added  during  the  modification  for  the 
type  2  design  approach  eliminated  air-entraining  vortices  and  was  partially 
responsible  for  the  headwater  elevation  lowering  from  1203.6  to  1202.5.  The 
1 0 . 5 - f t - radius  belimouth  first  tested  on  the  type  6  design  approach  was  par¬ 
ticularly  effective  in  preventing  debris  accumulation  since  the  radius  allowed 
the  debris  to  turn  or  pivot  into  the  tunnel.  The  use  of  a  bellmouthed  en¬ 
trance  is  therefore  recommended  whether  the  type  6  or  the  tvpe  7  design  piers 
a  r  e  u  s  e  d . 

29.  Of  the  seven  designs  studied,  the  type  6  and  type  7  design  ap¬ 
proaches  have  the  most  hydraulically  efficient  pier  configuration  for  passing 
debris  (Figure  8).  Both  types  resulted  in  maximum  headwater  elevations  of 


Figure  8.  Types  6  and  7  pier  dimensions 


approximately  1203.0  during  the  SPF  hydrograph  with  debris.  The  type  6  design 
had  some  debris  accumulation  at  the  entrance  to  Tunnel  C,  which  likely  was 
caused  bv  the  protruding  bridge  pier  on  tunnel  pier  2.  The  type  7  design  had 
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• debris  .uvuir.il.tt  ion  a;  the  cut  ranee  during  the  debris  hydrograph  test 

o'.  blunt  edges  and  flat  surfaces  at  or  be  Low  the  water  surface  tend 
;  a  .  a;;,  ;  u r ! ' '  i !  t  ■  i  ’.c  e  and  to  gather  debris.  Examples  of  this  were  observed  in 

t  :■  ;  pot  * a  I  tace  designs  in  types  1  through  b .  The  sloping  noses  with 

:  r.  •  >r  shaped  seel  ier.s ,  found  in  types  6  and  7,  helped  prevent  debris  from 
lodging  !>v  el  iminat  itig  a  "flat"  surface  for  the  debris  to  hang  upon  and  by 
logo 'log  v  aria  nts  that  tended  to  move  the  debris  toward  a  tunnel  opening. 

1 .  While  ttie  debris  tests  on  types  A  and  5  designs  did  not  result  in 
vert  orp 1 ng  <>t  the  diversion  structure,  these  designs  are  not  recommended 
oo.lt  -o,  v  n.otv  careful  debris  inventory  is  conducted  and  further  tests  are 

t.d  A  field  trip  was  taken  in  April  1988  to  inventory  debris  in  the 
lo.a:  Fork  floodplain.  Recommendations  pursuant  to  this  investigation 
ir.c’.ode  the  removal  <>f  all  housing  in  the  floodplain,  and  the  selective 
1.  i ng  of  a  1 !  trees  near  the  channel  bank  from  the  first  upstream  bridge 
o  wo  to  the  tunnel  entrance. 

i},  The  selection  of  the  design  approach  was  based  on  the  ability  of 
•  he  ootif i gurat ion  to  pass  the  debris  through  the  tunnels  and  on  downstream, 
i'i.t-  t  >;it  channel  should  therefore  be  assessed  for  its  capability  to  maintain 
the  project  flood  while  passing  this  debris.  A  bridge  located  approximately 
mile  downstream  of  the  exit  portal  is  a  potential  location  for  the 
aeeuirulat ion  of  this  debris.  Precautions  should  be  taken  to  prevent  this  from 
■:  >  v.-  c  u  r  r  i  n  g . 

31.  This  study  presents  a  topic  which  is  critical  to  the  safe  operation 
and  performance  of  many  hydraulic  structures,  that  is,  the  potential  hazards 
of  floating  debris.  While  research  studies  have  been  undertaken  to  assess 
structural  alternatives  such  as  booms  and  ’’debris”  basins  (Perham  1  987), 
further  research  is  still  ir-cilcd  in  debris  transport,  particularly  quantifying 
the  volume  as  it  relates  to  distribution  and  time. 
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Summary  Comparison  of  Data 
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Table  2 


Debris  Distribution 


1  tern 

Desc  r iDt ion 

Model 
ft  ^ 

Prototype 

ft.3 

Large  Items 

Large  trees 

i'D  to  90  ft  tall 
(includes  28 
fabricated 
trees  ) 

14 . 3 

386,100 

Subtotal  (20%  total  volume) 

14.3 

386, 100 

Small  Items 

Small  trees 

20  to  35  ft  tall 
( includes 
variable- 
diameter 
logs  30  ft 
long) 

20.0 

340.000 

Branches/twigs 

5  co  20  ft  tall 

20.0 

540.000 

House/trailer 

Siding 

10  bv  10  ft 

15  bv  10  ft 
( variable 
thickness 
and  buovancv) 

5.5 

12.0 

148.500 

324.000 

Subtotal  (80%  total  volume) 

57.5 

1,552. 500 

Total 

71.8 

1.938.600* 

*  Total  volume  is  less  than  1  percent  greater  than  the  estimated 
volume  of  1.92  5.000  1 1 ^ . 


Table  3 

Entrance  Losses 
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2.19 

Type  7 

wo/debris 

1201.9 
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0.01 

2.19 
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